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The thick brane model is considered in f(R) ~ R n gravity. It is shown that regular asymptotically 
anti-de Sitter solutions exist in some range of values of the parameter n. A peculiar feature of this 
model is the existence of a fixed point in the phase plane where all solutions start, and the brane 
can be placed at this point. The presence of the fixed point allows to avoid fine tuning of the model 
parameters to obtain thick brane solutions. 
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■ I. INTRODUCTION 

s : 

At the present time, the study of the structure and evolution of the Universe is at an interesting stage. The 
, necessity of a consistent description of the present Universe demands the creation of a unified theory of elementary 
particles and cosmology. A very promising way is the consideration of models of the Universe in higher-dimensional 
theories. Such investigations were initiated in the works of Kaluza and Klein in the 1920s for the unification of the 
O" 1 two fundamental interactions known at that time - gravitation and electromagnetism - within the framework of a 
5^ ' unified five-dimensional theory. Later on, similar ideas were used for the unified description of the four currently 
known fundamental interactions within the framework of superstring theories with several extra space dimensions. As 
in the case of the Kaluza-Klein theories, in superstring theories, it is supposed that our four-dimensional space-time 
^vq , results from the spontaneous compactification of a higher-dimensional space. 

^ ' At the same time, models of the Universe with non-compact (and even infinite) extra dimensions are under consid- 
er | eration [l[ (for a review, see also @, Q). In such a theory, it is supposed that we live on a thin leaf (brane) embedded 
into some higher-dimensional space (bulk), and matter is somehow confined (trapped) on the brane. The existence 
of extra dimensions then allows to resolve a number of old problems in high-energy physics (such as the problem of 
. . mass hierarchy, stability of the proton, etc.). 

All branes can be divided into two classes: thin branes and thick branes. In the first case, one has a delta-like 
localization of matter on the brane [l|. From a realistic point of view, however, the brane should have some thickness. 
The inclusion of a brane thickness then yields new possibilities and new problems (for a review, see, e.g., Q). Such 
a brane must satisfy two major requirements: 1) the solutions should be regular and asymptotically flat, or de Sitter 
ones (or anti- de Sitter ones); 2) ordinary matter should be confined to the brane. 

Most thick brane models employ scalar fields within the framework of Einstein's theory of gravity (see e.g. the 
review 4] and references therein). However, one might expect the existence of brane- like solutions also for some kinds 
of modificated gravity theories, the so-called higher-order gravity theories (HOGT). In such theories the action of the 
Einstein-Hilbert gravitational Lagrangian is supplemented by further terms, which are curvature invariants [5| . (Such 
a modification is based on the effect of the interaction of quantum matter fields with the classical gravitational field.) 
This allowed to avoid an initial cosmological singularity and to construct regular cosmological models of the early 
Universe [6| . Later it was shown that in such type of models a stage of inflation can exist [7| . 

Currently, this last possibility is widely used for the description of the present accelerated expansion of the Universe. 
This acceleration can be explained by the presence of some antigravitating substance - the so-called dark energy. The 
description of dark energy can also be realized within the framework of f(R) theory, where f(R) is some arbitrary 
function of the scalar curvature R. By choosing f(R) ~ R n , it was shown that such models are in good agreement 
with several different sets of observations On the other hand, such theories can be successfully employed for 

the description of dark matter as well. 
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HOGT with more complicated combinations of curvature invariants are also under consideration. In particular, in 
the low-energy limit of M-theory the Gauss-Bonnet invariant appears 

G = R 2 — iR^R^ + R^paR tll/pa ■ 

It was shown that such models, on the one hand, do not contradict observations within the Solar System and, on the 
other hand, successfully describe the present accelerated expansion of the Universe [14| . These models can be used 
in the description of an effective equation of state both for an effective cosmological constant and for the dynamic 
case (quintessence, phantom dark energy), and also for the description of the transition of one type of dark energy 
(quintessence) into another one (phantom energy). Also there are theories which employ both f(R) and Gauss-Bonnet 
terms to describe dark energy [la ]. 

Another usage of HOGT consists in the consideration of higher-dimensional cosmological and astrophysical models. 
In particular, in Ref. [lH brane world and black hole models with higher-dimensional action 

aR 2 + bR^B!" + cR^peR^P" + ^R-A + L m , 

(where L m is the matter Lagrangian; a, b, c are arbitrary constants) were considered. The obtained results allow to 
estimate general properties of models within the framework of HOGT. 

Furthermore there are some results obtained by the application of HOGT for the creation of brane-world models 
In particular, Parry et al. [l8[ considered a brane model in f(R) = R + aR 2 theory. Making use of the conformal 
equivalence of such gravity models and Einstein-Hilbert gravity with a scalar field, the authors rewrote the f(R)- 
equations in the form of the Einstein equations with some scalar field source. They showed that in such models 
brane-like solutions exist. 

Whereas Parry et al. [l8[ considered thin branes, in this paper we investigate thick branes in such f(R) theories to 
see, whether this leads to new and physically more acceptable results. 




II. EQUATIONS AND SOLUTIONS IN f(R) ~ R n THEORY 



We will work in a five-dimensional spacetime. The corresponding gravitational action can be taken in the form 

"!+/(*)] - W 

where f(R) is an arbitrary function of the scalar curvature R. (Here we employ units such that 8t:G — c = I.) 
Variation of the action ([T]) with respect to the 5-dimensional metric tensor g A B led to gravitational equations: 

R B A -\s B A R = n, (2) 
where capital Latin indices run over A, B, ... = 0, 1, 2, 3, 5, and 

*° = " { (1) *t- 5* + - ^ M > {%) J <*> 

specifies the effective geometric matter source with the nontrivial dependence on curvature; the semicolon denotes 
the covariant derivative. One can see that the gravitational equations in /(i?)-gravity rewritten in the form ^ have 
a structure that coincides with the standard general relativity equations when the source of the gravitational field is 
the effective energy-momentum tensor ([3]) . One can check that the energy- momentum conservation law is satisfied as 
well (see, e.g., [gj]). 

Here we will focus on a special choice of f(R) in the form 

f(R) = -aR n , (4) 

where a > and n are constants. As shown in Refs. Sj, where the present accelerated expansion of the Universe 
was considered, there are some ranges of n which do not contradict the observational cosmological data. Therefore it 
seems natural to consider these values of n for brane models as well. 
We will adopt the flat brane model with the metric 




ds 2 = e 2y ^r] al3 dx a dx - dz 2 , 



(5) 
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where the warp factor function depends on the fifth coordinate z only, and r) a p — {1,-1,-1,-1} is a Minkowski 
metric. Inserting this metric into Eqs. @ and (J3J, one obtains from the (*) component of the Einstein equations 



cPp 
1 dy 1 



dp 
op— — 



1 



dy 32p 2 f RR 



4p 



dp 
dy 



- P ) fn - ^/ ~ 6 P 2 



(6) 



where the new variable p — dy/dz was introduced, and the index R denotes the derivative with respect to the scalar 
curvature R. Eq. ((6]) is a third-order differential equation with respect to the metric function y, and the remaining 
components of the Einstein equations are fourth-order ones. Using the expression for f(R) from ((4]) one can obtain 
the equation for y in the form 
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n y' 
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(7) 



where the prime denotes the derivative with respect to z. We note, that by introducing the scaled variable z = az 
1 

with a — a 2 ( 1 -™) , Eq. (J7|) becomes independent of a. Thus it is sufficient to solve the equation for a = 1. All other 
solutions are obtained from this solution by scaling. 

One can also see from Eq. ([7]) that the first derivative y' cannot take the value zero unless y"(z) — as well. As will 
be shown below there is only one point in the phase plane where both y' and y" are equal to zero - the fixed point. 

Due to its complexity, we will seek numerical solutions of Eq. ([7]). But before, let us investigate the qualitative 
behavior of solutions of Eq. ((7|). For this purpose, we rewrite it as a system of three first-order differential equations 
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= p, 
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= v, 
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n p 





(8) 
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The fixed point of the system is 



2 n(n - 1) a8 n n(n - 1) 



A = {p^ Q,v -> 0}. 



v+-p 



(9) 



In order to analyze the behavior of solutions at the fixed point, let us consider a variation S -C y in the neighborhood 
of the fixed point for Eq. (J?) . The corresponding equation for 5 then becomes 



n y' 



6" 



n y n 



1/3 



-5' = 0. 



(10) 



Taking into account that at the fixed point denoted by z = Zf p the derivatives vanish, i.e., y'^ p = and y'j p = 0, we 
seek a solution of Eq. ([7]) in the neighborhood of the fixed point of the form 



y = yf P + y'l'p 



(z - z Sp f 



Then Eq. JT0| takes the form 



with the solution 



■5" + 



n{z - Zfp) n{z - z fp ) 



S' = 



S = S a (z-Zfp) 2+4 ^. 

To be sure that this solution decays faster than y, it is necessary to suppose that 

4 

2+->3 0<n<4, n^l. 

n 
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Let us next analyze the behavior of solutions in the plane {v,p} near the fixed point following the approach of [l| 
As an example, let us consider the important case with n = 2 (and a — I). In this case Eq. (J7J) takes the form 



p" + f(p,p')p' + g(p) = o, 
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or, equivalently, 



where 



P = v, v' = -g(jp) - f(p, v)v, 



f(P,p') = -~ + 4p, 5O) = ^p 3 - ^p. 



It is obvious from the above expression that 



±00 



J g(p)dp = G(± 



00 = +00. 



This is a necessary condition to be satisfied in the analysis of such systems [19] . According to [1J|, in this case the 
origin of the coordinates is then a single fixed point. It is surrounded by the so-called "curves of energy" with the 
equation 

w(p, v) = GO) + \v 2 = C. 

All of these are closed curves. In our case the function g(p) is odd, and therefore these curves are symmetrical with 
respect to both coordinate axes. 




FIG. 1: The phase portrait for the case n — 2 (left panel) and n = 4/3 (right panel), a = 1. A is a repulsive node, B, C denote 
the asymptotic points (|12|) at z — > T°°> respectively. The different curves correspond to different initial values of the variable 

y- 

The behavior of the solutions in the neighborhood of the fixed point (and the curves of energy) can be estimated 
as follows: from the second and third equations of the system (JHJ) , it is possible to obtain an equation near the fixed 
point 



dv 
dp 



1 v 
2p 



where D is an integration constant. This is a set of curves escaping from the node p = 0,v — 0. Specifying the 
boundary conditions in neighborhood of the fixed point A., one can then draw a phase portrait of the system ([8]). 
presented in Fig. [TJ As a further example, the model with n — 4/3 is also presented in this figure. 
The asymptotical form of the solution for arbitrary n is given by 



£/oc 1 ^ j 3 



12 



ail 



(12) 



One can see that there is an upper bound for the parameter n, n < 5/2. Taking into account the condition from (jTTJ) 
one can conclude that such a type of thick brane solution can exist only in the range < n < 5/2, n^l. 

Let us now discuss the corresponding distribution of the effective energy density Tj = —3 (y" + 2y' 2 ) for the case 
n = 2, shown in Fig. [5] Asymptotically Tq goes to a constant negative value: ^(±00) ~^ ~^v'oo = — 0.9, and the 
5-dimensional scalar curvature — 8y" + 20y a goes to a positive constant, ^/±oo) = ^' Thus this corresponds to 
an asymptotical anti-de Sitter solution. 
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FIG. 2: The effective energy density Tq for the case n = 2, a = 1. The point Tq = has been placed at z = by shifting 
z — ^ « — Zo. 



III. TRAPPING OF MATTER 



In this section we consider the trapping of a test scalar field on the brane considered above. For this purpose, let 
us use the approach suggested in the paper [2(| ■ We consider the test complex scalar field \ with the Lagrangian 

t 1 a * 1 2 * 

L x = o^aX 9 x ~ 2 m ° X X ' 
where mo is the mass of the test field. Using this Lagrangian, we find the equation for the scalar field 

7^ (,^ sAB ^) " (13) 

Here x is a function of all coordinates, x — xi xA )- Taking into account that the canonically conjugate momenta 
= (E,p) are integrals of motion, we will seek a solution in the form 

x (x A ) = X(z)exp(-ip li x^). 
Inserting this ansatz into Eq. (fT5)l . leads to the equation for X(z) 



X" + v/^WPm - ml)X = 0, 

or, taking into account that = e~ 2y (E 2 — p 2 ), we obtain 

X" + [(E 2 - p 2 ) e 2y - m 2 e^] X = 0. 

According to Eq. (|12|) . asymptotically = k n \z\, k n > 0. That is why the dominant term in the above equation 
will be the term with e v . Thus 



X" - m 2 e ik ^X = 



with the asymptotically decaying solution 



where C is an integration constant. 

As a necessary condition for the trapping of matter on the brane, one should require finiteness of the field energy 
per unit 3-volume of the brane, i.e., 

oo oo 

E tot [ X }= J T*^~gdz= J e Ak ^ [e~ 2k ^ (E 2 + p 2 )X 2 + m 2 a X 2 + X' 2 } < oo, (15) 



6 



and also finiteness of the norm of the field x 



oo 




■OO 



e U n \z\ x 2 dz _ 



— oo 



— oo 



From the solution (|14l) it is evident that both E to t and ||x|| converge asymptotically. 



IV. CONCLUSION 



We have considered the 5-dimensional thick brane model in f(R) ~ R n theory. It was shown that regular solutions 
with asymptotic (II 2[) exist in the range < n < 5/2, n ^ 1. Space is asymptotically anti-de Sitter in such a brane 
model. 

A special feature of the model is the presence of the fixed point A. Since Eq. ([7]) is invariant under the shift of the 
independent variable z z + zq, the position of the brane is arbitrary and it can be placed at any point on the axis 
z, including z = 0. The most natural assumption seems to be that the brane is situated at the fixed point A at z = 0. 
Then y(0) = const, y'(0) = 0, y"(0) = holds on the brane. 

The presence of the repulsive fixed point A (node) in the system allows to set the boundary conditions arbitrarily 
in the neighborhood of the fixed point. Anyway, all solutions will escape from A and tend to the asymptotic value 
(|12[) . It allows to not provide any special conditions for the model parameters (fine-tuning conditions) that is typically 
necessary for other brane models (see, for example, Ref. (20|V 

Consideration of the behavior of a test scalar field in the bulk has shown that such a field is trapped by the 
/(i?)-brane. Note that the trapping is purely gravitational. 

For n = 2, this f(R) theory has also been employed for the thin brane model [Hj]. In the vicinity of the brane 
(located at z — 0) a similar behaviour of the metric function is found, y(z) ~ z 3 . However, at some finite value z s , 
a singularity is encountered. Thus the thin brane model appears to be plagued by bulk singularities, very much in 
contrast to the thick brane model considered in this paper. 
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